Fast and robust quantum computation with ionic Wigner crystals 
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We present a detailed analysis of the modulated-carrier quantum phase gate implemented with 
Wigner crystals of ions confined in Penning traps. We elaborate on a recent scheme, proposed by 
two of the authors, to engineer two-body interactions between ions in such crystals. We analyze for 
the first time the situation in which the cyclotron (uj c ) and the crystal rotation (oj r ) frequencies do 
not fulfill the condition ui c — 2cj r . It is shown that even in the presence of the magnetic field in 
the rotating frame the many-body (classical) Hamiltonian describing small oscillations from the ion 
equilibrium positions can be recast in canonical form. As a consequence, we are able to demonstrate 
that fast and robust two-qubit gates are achievable within the current experimental limitations. 
Moreover, we describe a realization of the state-dependent sign-changing dipole forces needed to 
realize the investigated quantum computing scheme. 

PACS numbers: 03.67.Lx,37.10.Ty,37.10.De,45.50.-j 
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I. INTRODUCTION 

Despite the huge experimental progress to cool, trap, 
and manipulate single particles such as atoms and 
molecules at the quantum level, the way to build up a 
quantum computing hardware working with several hun- 
dreds of quantum bits (qubits) in a coherent and control- 
lable manner is still long. By means of quantum optimal 
control techniques it is possible, at least theoretically, to 
perform one- and two-qubit quantum gates with fidelities 
above the demandingthresholds of fault-tolerant quan- 
tum computation [iMsj]- These thresholds fix an error 
between 0.01% to fractions of a percent [Io|, Hlj . Up to 
now, only with cold trapped ions quantum gates with a 
fidelity of 99.3% have been experimentally demonstrated 
EH [l3| . which is not too far from the aforementioned 
thresholds. Similar fidelities have been also obtained for 
small quantum algorithms [l4l . Il5j . 

Nowadays, however, most of the experimental efforts of 
the atomic and molecular physics community are concen- 
trated in the design and fabrication of microtraps, both 
for ions [H, [I?} and neutral atoms [HI, [lji| . Even though 
these efforts are important, significant technical issues 
related to the miniaturization and trapping methodolo- 
gies arise when scaling to many particles, and therefore 
new strategies have to be devised. A possible solution 
to the problem is the separation between the qubits used 
as quantum memory and the ones employed to process 
the information [2(| or, alternatively, the exploitation of 
quantum distributed networks [2l| . Another approach, 
instead, consists in the use of collective states of atomic 
ensembles with a multilevel internal structure as qubits 
0. 

Apart from these technological efforts and alternative 
solutions, nobody can yet say which of the various phys- 



ical implementations will be the successful one. It is fair 
to say, however, that ions represent a good candidate 
to implement a multi-qubit quantum processor. Indeed, 
two-qubit gates with ions can be realized in about few 
tens of 11s [23|, and qubits stored in internal elec- 
tronic degrees of freedom of an ion have coherence life- 
times ranging from 1 s to 100 s or more [24j . 

Coulomb — also named classical Wigner — crystals 
confined in Penning traps are natural candidates for 
a quantum memory, since the separation among ions, 
about 10 /zm, allows to individually manipulate their in- 
ternal degrees of freedom. Such a trap scheme uses static 
electric fields to confine charge particles in the axial direc- 
tion (the z axis in Fig. Q]), whereas the radial confinement 
is provided by a strong uniform magnetic field along the 
axial direction. Currently, Penning traps allow to trap up 
to 10 8 ions [25[. An appropriate choice of the trap pa- 
rameters (e.g., tight axial confinement) allows the ionic 
ensemble to crystallize in a two-dimensional (2D) hexag- 
onal lattice configuration with an inter-particle spacing 
on the order of tens of /xm [26| , and therefore to manipu- 
late a large number of qubits without specific micro-trap 
designs. The high phonon mode density, however, does 
not permit to resolve single modes for sideband cool- 
ing. Hence, Doppler and sympathetic cooling are the 
most natural techniques to be employed; we also note 
that Sisyphus cooling might be an alternative methodol- 
ogy [13 ■ Current experiments, however, performed with 
Doppler cooling, can reach temperatures of few mK [28j , 
that is, a high thermal occupation number distribution 
of phonon modes. Nonetheless, efficient quantum com- 
putation and production of small cluster states are the- 
oretically possible [IH [3(j, and recently full control of 
the qubit Bloch vector with ^99.85 % fidelity for Rabi 
flopping has been experimentally demonstrated |3l| . 

The two-qubit gate scheme considered in the propos- 
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FIG. 1. (Color online). Two-dimensional Coulomb crystals 
of ions in a Penning trap rotating at frequency oj r . To manip- 
ulate the internal states of the ions, laser beams can address 
single sites or multiple ions. 



als of Refs. [29|, |30( is based on the so called "pushing 
gate" (or its variant, the modulated-carrier gate), where a 
spatially inhomogcncous laser field together with an ap- 
propriate combination of polarizations and frequencies 
induces a state-dependent dipole force on two nearest 
neighbours of the 2D Coulomb crystal (see Fig. [TJ. De- 
pending on the configuration of lasers and polarizations 
the displacements of the ions away from their equilibrium 
positions can be either perpendicular to the plane of the 
crystal 1231 or along the in-plane separation of the ions 
[3fl [32H34| . The coupling between these displacements, 
mediated by phonons, yields entanglement of the internal 
states (qubits) of the ions, that is, the desired quantum 
gate between ions. 

In addition to the confinement, the radial electric and 
axial magnetic fields induce a drift that causes in-plane 
rotation of the crystal (see Fig. [1]), whose frequ ency 
w r /(2-7r) is typically on the order of few tens of kHz j26| . 
There are two possible solutions to our quantum hard- 
ware design: either we use a co-rotating (with the crys- 
tal) laser beam in order to realize the desired two-qubit 
quantum gate, or we have to perform the gate in a time, 
r g , such that the crystal rotation has a negligible effect 
on the gate operation. The latter solution translates in 
the condition w r r g /(27r) <C 1. Such a requirement is in- 
strumental because, in order to accumulate the necessary 
two- ion phase for the quantum gate we aim to implement, 
the ions have to experience the applied light force for the 
entire gate operation, or else, the required phase would 
be achieved only partially. 

While the former solution applies for all rotation fre- 
quencies, but relies on a more sophisticated experimental 
setup, the latter restricts the range of possible values of 
ui T . Thus, in both proposals [29|,[30|, where the rotation 
and cyclotron frequencies fulfill 2 uj t = uj c , the aforemen- 
tioned condition is satisfied when r g is on the order of 
ns, whereas the modulated-carrier gate of Ref. [3(| had 
T g = 5 ijls. Given the above, such a proposal requires a 



co-rotating laser beam. Thus, by maintaining 2 ui x = uj c 
one should reduce lo c . This approach, however, would 
not help since the smaller the cyclotron frequency is, the 
longer the gate operation. Instead, if we abandon the 
assumption 2 tJ r = cj c and look at moderate rotation fre- 
quencies, at the expenses of possible large modulations 
of the force, we are able to fulfill Ld r T g /(2ir) <C 1. Addi- 
tionally, low rotating frequencies result in low densities 
and large inter-particle spacing, and therefore in an easier 
way to address the trapped ions with a laser field. 

Thus, the main goal of this work is to analyze this 
regime and, at the same time, to perform robust two- 
qubit gates within a range of experimentally achievable 
temperatures. 

In the following we shall we present the general theory 
of the modulated-carrier push two-qubit gate (Sec. HIT) 
with details that were briefly mentioned in Ref. [301 ] . 
Subsequently, in Sec. IIII1 we investigate the situation in 
the presence of the magnetic field in the rotating frame of 
reference and the relative gate performance. In section lTVl 
we describe how to physically realize the state-dependent 
force required for the proposed quantum processor, and 
Sec. fVl summarizes our results and provides some future 
prospectives. 



II. MODULATED-CARRIER GATE 

In the following we make the approximation that the 
Wigncr crystal is a rigid body, which is a good approxi- 
mation in the magnctohydrodynamic regime (one compo- 
nent plasma) or at equilibrium (35j . Hence, in a rotating 
frame, the Hamiltonian of a crystal with N ions written 
in cylindrical coordinates [r = (r, 8, z)] is given by (35j 
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with m being the mass of the ion and 
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(2) 



Here T is a parameter describing the trap geometry and 
applied voltage on the electrodes [36] , eo the vacuum per- 
mittivity, e the electron charge, and u c = eB/m is the 
cyclotron frequency. We see, from the first line of Eq. ((T|), 
that there exists a special rotating frame, u = u> c /2, such 
that the minimal coupling disappears, and, in this sec- 
tion, we shall consider such a frame of reference together 
with uj t = w c /2 (i.e., the frame of reference coincides with 
the crystal). We note, that with "minimal coupling" we 
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refer to the interaction p-A. Such terminology is typically 
used in quantum field theory [37j ■ 

Finally, it is worth to remind that the gate we aim to 
accomplish realizes the true table |ei, e 2 ) — > e tfleie2 |ei , €2) 
with £1,2 = 0, 1 and = 6 m - 9 01 - 9 10 + 9 n Hi |Hj|. 
Specifically, we arc interested in a phase gate with 9 = tt, 
which, up to additional singlc-qubit rotations, is tanta- 
mount to a two-qubit controlled NOT gate [24j . 



A. Normal modes and canonical quantization 

The Hamiltonian ([TJ in cartesian coordinates [r 
(x, y, z)] reduces to 



N 

=£ 

k=l 



2m 



,2 _2 



v c , 

(3) 



where uj z = y 4T/m is the axial frequency, and u! xy = 
1/2(^2 - 2^2)1/2 t h e in-plane one. 

By performing a Taylor expansion of the potential up 
to second order around the stable equilibrium configu- 
ration, obtained by minimizing the total crystal energy, 
we can express the Hamiltonian in the new coordinates 
Qn.r) = T) n — ij®, that is, the displacements from the equi- 
librium positions. Hence, it is possible to determine an 
orthogonal transformation M such thalQ 



Hr 



p2 

2m 



m 2 2 



(4) 



with Q n>ri = J2k,n M n,n;k,nqk,iJ.[= Mq], and P = p. 

Now, we perform the canonical quantization and we 
introduce the creation (annihilation) operators a K (clk) 
for each mode K = (n,?7), along with the harmonic os- 
cillator ground state size uk = \/ h/muJK- Hence, the 
(phononic) Hamiltonian operator reads 



Hr = ^2 frw K (a\ ( a K + 1/2), 



(5) 



K 



where for the sake of simplicity wc drop y^f) in Hr. 

B. Adiabatic and oscillatory quantum gates 

Let us consider a spatially inhomogeneous laser field 
appropriately detuned from the internal states such that 



it produces a state-dependent displacement of the ions. 
Then, the matter-field interaction, in the electric dipolc 
approximation, becomes 



N 



K 



V2" 



(0) 



where fj is the three dimensional force due to the gradi- 
ent in the laser intensity, and crj is the z Pauli matrix. 
Here the following relation for the displacement coordi- 
nate operator 



qi< = J2 M k>-,kQk> = J2 MK>-,K^=(a K + a K ) (7) 

K' K> ^2 

has been used. Thus, we have [K = (j, /x)] 



(8) 



where [fj'it)]^ is the \x' = x,y,z component of the three- 
dimensional vector fj'{t). Hence, the full problem re- 
duces to 3iV independent, driven oscillators. 

When the temporal profile of the force fulfills the 
condition lim t _ i .± 00 /(t) = 0, the unitary time evolu- 
tion operator is given by Uxit) — e~ l $ K & cxp(f3Ka K — 
PVclk) exp(— iuixt a^aft-), where cf>K and /3k satisfy the 
differential equations [38|, l39j | 



(XK 

p K = -iw K p K +ij-j=f K {t), (j) K 



hV2 



f K (t)Re[f3 K (t)}. 

(9) 

Given that, let us consider the adiabatic regime regime 
where fx it) varies slowly with respect to ujk [H}. Adi- 
abatic elimination, by taking (3k ~^ 0, yields 



UK.fK(t) 

HukV^ 



5 A" 



it) 



^W-ujk 



(10) 



Thus, the displacement of a normal mode K induced by 
the gate is proportional to the force applied, and can 
be made zero independent of the initial phonon state 
by starting and ending with zero force. This elimi- 
nates any potential error due to entanglement between 
phonons and the internal states of the ions. Similarly, 
the overall phase accumulated ^ K 4> K i T e) does not de- 
pend on the initial phonon state. However, for a gate 
occurring over a time interval [0, r g ], the final qubit state 
has applied exp(— i J2 n j '^njO'n^'j)^ where the two-particle 
phases arise from 



1 Hereafter we shall use latin symbols for index the ions and greek 
symbols for the cartesian coordinate of the force vector acting 
on the ions. 



/£(*)= E M K .^MK-, n , n [fM*l?n(t)]r,%K- (11) 
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Thus, the two-particle phase is given by 



! "j 



J2 S m j) HdtimUfrMn, (12) 



where the term outside the integral is a shape indepen- 
dent form factor, whose specific form is given by 



c("i) 

[IT] 



E 

A" 
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K 



2h 2 



■M, 



10 K 



K;j,iiMK;n,T)- 



(13) 



Hence, we can think about (|12p as a convolution of the 
forces on the two particles, modified by the form fac- 
tor representative of the characteristic oscillator variance 
over its frequency, which is overall proportional to w^- 1 . 

Now, let us consider a scheme with a force f(t) —> 
cos(vt) f (t) , where the carrier frequency v must be much 
larger than the modes of frequency ujk that are coupled 
to the force (this averages out any net displacement). If 
the modulation f(t) is slow as compared to v (but with no 
restriction with respect to ujk)-, we can perform a similar 
adiabatic elimination as above, and get a gate with the 
same desirable properties that can operate non-trivially 
on arbitrarily in-plane vibrational modes at very high 
temperatures. 

For adiabatic elimination with respect to we choose 



the Ansatz (3 K = (3 



K ' 



for each mode. By 



inserting this Ansatz into the differential equation © we 
obtain 



+ i 



2 V / 2V 

0^f K (t)-i(uJ K + v)P+. 



(14) 



Separate adiabatic elimination of (3 K and f3 K yields /3 K = 

&KfK(t)/{2V2h(LJK ± v)]- As before, in the pure adia- 
batic regime, we find that the displacement of a normal 
mode induced by the gate is proportional to the force ap- 
plied. Again, it can be made zero independent of the ini- 
tial phonon state by starting and ending with zero force, 
and therefore eliminating any potential error due to en- 
tanglement between phonons and the internal states of 
the ions. 

Now, we examine the two-particle phase induced in 
this new scenario. The time evolution of the phase is 
governed by 



'AT 



K 



2h 2 (o 



K 



V 2 ) 



cos 2 (vt)f K {t), 



(15) 



where the quickly varying component cos 2 (vt) can be re- 
placed with 1/2. As described in the adiabatic regime, 
the overall phase accumulated (^K^g), for a gate oc- 
curring over a time interval [0,r g ], does not depend on 



the phonon initial state. In this case the pulse-shape 
independent form factor is given by |3Q] 



a 2 K uj K 



K 



\h 2 {v 2 



-M 



K;j,iJ.MK;n,ri 



(16) 



J K1 



Performing a Taylor expansion in 1/v 2 the first term is 
proportional to J2k K j ,, ^1 k ,, ,, = <5j>£|»,»j {5j,n indi- 
cates the Kronecker symbol). This follows from the fact 
that M is an orthogonal matrix. Physically, this arises 
due to the coherent averaging of in-phase oscillating ions 
without any virtual excitation of phonons — accordingly, 
no two-body phase should be expected. The second term 
of the expansion is non-zero and yields 



UkMk^Mkm + O (V 



(17) 



K 



Compared to adiabatic push gates, the modulatcd-carrier 
gate is inverted in sign and it is multiplied (in phase) by a 
factor (lok/v) a /2 [see Eq. ([T3j)], In the case of a lateral 
operating modulated-carrier gate with oj xy <C v <C uj z , 
the accumulated phase is enhanced by a factor (lo z /v) a /2 
with respect to an adiabatic push gate with a force mov- 
ing the ions in the axial (z) direction for the same laser 
parameters. Given that, the gate time needed to per- 
form a 7r-phase gate is reduced. In the opposite case, 
that is, for an adiabatic in-plane push gate (u>k ~ U) X y), 
and for the same laser parameters, the lateral modulated- 
carrier gate is reduced in phase, and therefore a longer 
tv is required. Thus, compared to the proposal of Rcf. 
[291 ] , where the push gate operates in the axial direction, 
our modulatcd-carrier gate working with in-plane modes 
yields a larger two-ion phase for a given set of laser pa- 
rameters, and therefore it enables to perform a larger 
number of quantum gates within the coherence time of 
the system. 



III. MODULATED-CARRIER GATE WITH 
MINIMAL COUPLING 

In this section we analyze the situation where uj y =/= 
cj c /2, for which we have three reasonable choices for the 
rotating frame of reference: 

• F\ coincides with the lab frame, where the equilib- 
rium positions of the ions in the crystal are time- 
dependent and the minimal coupling does not van- 
ish; 

• Fi rotates with a frequency uj = w c /2, as in the pre- 
vious section, where the minimal coupling vanishes, 
but the equilibrium positions are time-dependent; 

• .F3 rotates with a frequency uj — uj r , where equilib- 
rium positions arc time- independent, but the min- 
imal coupling does not vanish. 
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A. Equilibrium configuration of the crystal 



B. Quadratic expansion of the Hamiltonian 



Let us discuss which of the frames of reference f 1,2,3 is 
more suitable to numerically determine the equilibrium 
configuration of the system for a fixed (a priori) value of 
total canonical angular momentum Pg0. Since we are not 
concerned with relativistic velocities, the electromagnetic 
fields involved in the problem are the same in all frames 
of reference. Consequently, the angular momentum of 
an ion in a frame rotating with uniform angular velocity 
with respect to the (inertia!) laboratory frame coincides 
with the one in the latter [40| ■ This conclusion allow us 
to find the equilibrium configuration of the crystal, for 
a given value of Pg, by choosing a frame of reference ro- 
tating with angular velocity u> = lo c /2 (the frame Fi in 
the above outlined list) in such a way that the coordinate 
systems at the initial time t = of F2 and F3 do coincide. 
Such a choice simplifies the numerical minimization pro- 
cedure, because the minimal coupling in the (classical) 
Hamiltonian vanishes. We underscore, however, that F2 
is utilized only at time t — for the determination of 
the equilibrium configuration of the crystal. Instead, for 
times t > we use F3, where the equilibrium positions 
are time-independent. With such a choice the numerical 
effort in order to assess the gate performance is signifi- 
cantly reduced. 

Besides this, we also note that not all rotation frequen- 
cies w r of the crystal allow to have a stable configuration, 
that is, ions confined within a well-defined spatial region. 
Indeed, by rewriting the addend of the second sum in 
Eq. (HJ as 



T(2z, 



w r )r 



<4 + Pri) 



(18) 



we see that the potential is confining if and only if j3 is 
positive. Here the anisotropy parameter /? is defined as 



a 2 z uj c 



1 - 



1 

2' 



(19) 



where a z = uj z /lu c . Importantly, /3 relies only on a z and 
the ratio aj r /uj c . Thus the range of admissible frequen- 
cies is: uj m < Lu r < uj c — w m , where w m = uj c j2 — u xy is 
the magneton frequency (35j . Of course, the admissible 
regime is also constrained by the condition a z < 1/ \[2. 
In order to access lower rotation frequencies, the trap pa- 
rameters might be changed by increasing u> xy , that is, by 
lowering lu z . Attention has to be paid, however, when 
uj c and uj z are changed, since due to such a manipulation 
different structural phase transitions may occur. In par- 
ticular, we are interested in the limit (3 <g 1, where a 2D 
hexagonal lattice structure appears [35j. 



2 When u) t ^ Wc/2, the total canonical angular momentum Pg ^ 0, 
but it is still a constant of motion 13511 . 



Let us introduce the typical scale of length £ s , momen- 
tum p s , and energy E s in our problem: 



\47reoTOWc 



p s = £ s muj c E s 



.(20) 



Then, the Hamiltonian (fTJ) in cartesian coordinates be- 
comes 



1 N 

h r {uj) = - 2J [pl k +p 2 Vk +pl k + {ykPx k - x k p Vk y, 



k=l 



1 N 

(l + 2a)] + z £ 



r.2 1 



1 N 

+ - lim y 

k,j= 



k=l 



k,j 



\ r k - rj + e| 



(21) 



where the substitutions H^(lj) — > Hh{uj)/E s1 {r,z) — > 
(r,z)/i s , (p x ,Py,Pz) (p x ,Py,Pz)/Ps, and a = oj/uj c 
have been introduced. The expression of the Coulomb 
potential, third line in Eq. (f2~Tj) . allows to obtain more 
compact formulae later in the present section. 

Next, given the equilibrium configuration (Fo,0) of 
each ion, we expand the Hamiltonian (|21j) to second or- 
der in the spatial displacement q = r ro and p around 
zero, namely 



H R (p, q) ~ H R (0, r ) + -d^d T , (22) 

where d T is the transpose of the row vector d = 
(qi,x,Pi,x,qi,y,Pi,v, ■ ■ ■ ,qN,z,PN,z), and H R = H R (0,r ) 
is the Hessian matrix. Its non-zero matrix elements are 
given by: 



d 2 H R 
dp 2 



1, 



d 2 H R 



d 2 H R 



dp Xk dyk dp Vk dx k 



2 H R 
drj k dfij 



= [1 - 2al + (6a 2 z - 1)6^] 

(1 - 6 kiS )[6 s j + (1 - ^,j)<J|fc- 3 -|,o] 
+ hm > ' ' — X 



x(-l) 



s=l 

<5fc,j 



\r k - r s + e\° 

(Vk - 1]s)(Pk - /is) 



\r k - r s + e| 
where r], fi = x,y, z, and k, j = 1, . . . , N. 
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C. Symplectic diagonalization and canonical 
quantization 

Hereafter we utilize the frame F3 that rotates at the 
frequency w r . Hence, we are allowed to drop H R (0, Tq) in 
Eq. (|22j) and the full Hamiltonian reduces to the 6N x 
6N- matrix H R {u) z ) = dH R d T /2. 

In order to perform the canonical quantization, we 
have first to transform the classical Hamiltonian H r (uj y ) 
in canonical form. A transformation S : (p,q) — > (P, Q) 
is canonical when the condition SIS T = J is satisfied, 
where J = i o v 0| ■ Since the Hessian matrix Ho 
is real and positive definite, Williamson's theorem [41| 
guarantees that 



SH R S T = W 



\ 



^3N 



(23) 



,3JV, 



where io k are real and positive numbers Vfc = 1, . 
and W is called the "Williamson form" of H R . 

Given that, we can recast the classical Hamiltonian as 



d; 



^ 3N 

— T£A* k>j a k + A k>j a 



(27) 



with A k j = S2k-i,j +iS2k,j, and where j is an odd integer 
[see the definition of the vector d after Eq. (|22|) ] . Then, 
the matter-field interaction has the following expression 



JV N-l 3 

V = ^2Wj ■ fj(t)]°j = + 1 ■ 77 3j+n(t)d2( n +3j)-l 

3=1 3=0 n=l 
3N 

= ^o£a fc + a k a{, (28) 
fe=i 

where J 7 = (fi, x ,fi, y , fi. z , /n.x, In, v , Jn.z), and 



Oik = 



^ 2V-1 3 
v 7=0 n=l 



(29) 



Thus the full Hamiltonian is: H = H R ((J T )+V = Y, k Hk, 
where H k = u k (a k a,k + |) + a* k a k + a k a\. 



j 3JV 3JV 
k=l fc=l 

where the new coordinates are determinated by the trans- 
formation A T = (S' _1 ) T d T . For the sake of simplic- 
ity, hereafter, we use the definitions Q k := A2k-i and 
P k := A2k Vfc = 1, . . . , 3iV. Thus, the Hamiltonian re- 
duces to 



H R {u t ) = -Y,"k{Ql+Pl), (25) 
fc=i 

that is, a sum of uncoupled harmonic oscillators. 

Similarly to Sec. Ill Al we perform the canonical quan- 
tization by promoting Qk,Pk to operators such that 
[Qk,P s ] = i8k,s- Besides this, we introduce the oper- 
ators a k = (Qk + iPk)/V2, a{ = (Qk - iPk)/V2 with 
[Sfc,aJ] = 5k, s- Hence, the quantized Hamiltonian is sim- 
ply given by 



H R (uJr) =Y,u k (aia k + -), (26) 
fc=i ^ ' 

and we note that the eigenvalues uj k are dimensionless. 

Finally, we rewrite the coupling between the ions and 
the inhomogeneous laser field. The displacement of the 
ion from its equilibrium position can be written as 



D. Two-qubit phase gate 

The time evolution of a phonon mode state, governed 
by the Hamiltonian H k , and a generic two-qubit state is 

|*k;*qbit(*)) = e-^WopfcWle-*^ |**;* qb it(0)) , 

(30) 

where £)[/3 k (t)] is the displacement operator [42j, @k(f) = 

-i /„' d S a fc (s)e Mfc ( s -*), and H° k = u k (a[a k + §) . 

In order to disentangle the external dynamics due to 
the phonons and the internal dynamics of the qubit states 
at the end of the gate operation, t = r g , the following 
condition has to be satisfied 39] 

X k = ^=r dte iuJkt a k (t) = Vfc. (31) 
V^fe Jo 

This condition, however, is more general than the adia- 
batic elimination we performed in Sec. Ill B[ whose aim 
was to highlight the difference in the accumulated two- 
particle phases among the most common quantum gate 
schemes based on pushing forces with off-resonant lasers. 

The necessary lateral force on the j-th and fc-th ion, 
\fj\ = \fk\ = A P hio xy cos(vt)e- t2 ^/\r° - f°\, is dc- 
tcrmined by setting the dimensionless parameter Ap to 
achieve a tt phase between the chosen pair of qubits. 
Then the fidelity is given by 
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F = min -j Tiy 



{Tr ph [(^ hit \iL(t) (p T (o) ® |$ qbit ) ($ qbit |){it(t) |$ qbit >] } 



= mm 
± 



n 



exp 



I _ e -Hu k /k B T 



,(32) 



where il(f) is the unitary evolution operator defined 
through Eq. (|3TJ|) . px(0) is the initial (canonical) density 
operator of the phonon modes at temperature T, |$ q bit) 
is the initial two-qubit state, and 



fa) fa) 



(33) 



ei,e 2 =0 



is the desired logical target state we aim to attain. The 
integral Iu q for q = 1 , 2 is given in Eq. (f3~T|) where the 
apex (j g ) refers to the ion we are considering, that is, 
j = j q in the sum of Eq. (|2"5|) . 

Since we aim to achieve r g -C 2-k/lo Vi we outline the 
following program: Firstly, we analyze the dependence of 
Lu r on the total angular momentum Pg. This is achieved 
by fixing a priori a value of Pg and then by determining 
the equilibrium configuration of the crystal, namely the 
positions and momenta of each ion (the most difficult 
part of the program) . Since the crystal is a rigid body, it 
holds p 9k = mrle k + er k A g (r k ) = mrl(uj r + oj c /2) [35[, 
and from this relation the rotation frequency is extracted. 
Such an analysis allows us to find the smallest value of 
u> T such that (3 < f3 c — 0.665/ VN is fulfilled, that is, a 
2D Wigner crystal configuration [35|. Then we choose 
the value of both w c and v in order to achieve high gate 
fidelity for a wide range of temperatures. 

The determination of the classical ground state is a 
multidimensional minimization constrained problem for 
which no deterministic and efficient algorithm is known. 
Here we used a variant of the Metropolis [43[ and the mul- 
tidimensional constrained Newton algorithm like the one 
of Ref. [44| ■ The first method allows us to sample ran- 
domly the relevant phase space region by choosing a slow 
decay of the acceptance probability and by using several 
annealing cycles. We then coarse-grained the obtained 
annealing trajectories into intervals, and we employed, 
for the lowest energy configuration on each interval, a 
Newton algorithm, which is very efficient in finding a lo- 
cal minimum provided that the initial value is already 
very close. We have checked the reliability of our numer- 
ical energy minimization for Pg = 0, that is lo t = 2lu c , 
by comparing the results of Ref. [44[ for the minimal 
excitation frequency for several numbers N of ions. 

We investigated the robustness of the modulatcd- 
carricr phase gate against temperature for a moderate 
number of ions N — 30 and a z = 0.70. In Fig. [2] the 
dependence of the crystal rotation frequency on the total 
canonical angular momentum is showed, whereas in Fig. 
[3]the gate infidelity for different values of the ratio r g /r r is 
displayed. The results of Fig. prefer to Pg = 4000 £jrmj c , 
for which we obtain the smallest value of |w r | in Fig. [21 
Beside this, we have for such a choice /3 = 3.4 x 10 ~ 4 , 



I 

whereas (3 C = 0.12, that is, a stable 2D hexagonal lattice 
configuration. Given that, Fig. [3] shows that in order to 
reduce by a factor 10 the ratio T g /r r the fast modulation 
frequency v of the force has to be (roughly) enhanced 
by a factor 10 as well. We also remark, that the three 
lines in Fig. [3] show an infidelity that is smaller for large 
gate operation times. The goal of the plot is to show 
how the modulation frequency increases when the ratio 
T g /r r is reduced for an infidelity smaller than 10~ 4 . Of 
course, by carefully tuning v one can easily get a smaller 
infidelity for faster gates. 

In the inset (left corner - top) the result of the gate 
infidelity for a cyclotron frequency 100 times higher is 
showed, that is, the same 7.608 MHz of the experiment 
of Ref. [28[ . Here there are two important features to be 
highlighted: firstly, the gate fidelity is more robust for a 
wide range of temperatures with respect to the previous 
case where ui c /(2ir) = 76.08 kHz has been considered. On 
the other hand, already for r g /r r = 0.1, the frequency v 
is on the order of hundreds of MHz. In the inset on the 
right (bottom) we show the gate infidelity again for the 
lo c =7.608 MHz but for a smaller modulation frequency 
v = 2.4 MHz that lies in the gap between the two bands of 
different radial modes (the so called E x B and cyclotron 
modes. See also Fig. [5]). In this scenario r g /r r = 10 
and therefore a co-rotating laser beam is required. In 
conclusion we see that for uj z ~ uj c / \[2 if we desire to 
avoid the employment of a co-rotating laser beam the 
only possible way is to achieve very high frequencies for 
the modulation of the state dependent force. 

Alternatively, one can consider a smaller value of a z , 
which basically shifts upwards the graph of Fig. [2J that 
is, by displacing the minimum of the |w r | closer to zero for 
large values of Pg. This is the situation depicted in Fig. @] 
for TV = 30, and a z = 0.02. Here it is possible to achieve 
gate operation times on the order of few fis with signif- 
icantly smaller values of the modulation frequency. In 
the figure v lies in the gap among axial and radial modes 
(see Fig. [5]). Furthermore with T g /T r = 6 x 10 -3 we do 
not need a co-rotating laser beam. This result is quite in- 
teresting since it works in a range of parameters that are 



currently employed in experiments (e.g., [31|). Finally we 



also note that in this scenario uj xy 3> u> z , which is oppo 
site to the requirement we identified in the case of larger 
a z when 2u> r = uj c . We note, however, that uj xv is not the 
actual radial frequency when uj c ^= 2uj t . Indeed, as shown 



in Eq. |T|), the centrifugal potential (i.e. 



,V2 



/2) 



modifies the confinement. Let us write u) T = (uj c — Soj)/2, 
where 8uj > 0. Then by substituting such definition into 
the second line of |T]) we obtain an effective radial fre- 



quency given by: ui[ 



eS 
xy 



parameters of Fig. S]we get w®*/(27r) 



- 2uj 2 z . With the 
31.47 kHz which 
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FIG. 2. (Color online). Ratio uj^/oJc vs. total angular mo- 
mentum Pg for a z = 0.7 and iV = 30. For Pg = we retrieve 
the well-known limit uj r — 2lo c , in which there is no magnetic 
field in the rotating frame. 



is significantly smaller than v z (~ 152 kHz), and there- 
fore the 2D lattice configuration is ensured. This fact is 
also confirmed by j3 = 4 x 1CP 2 < /3 C . 

Now, let us examine what is the required laser power 
in order to realize the gate and investigate the influence 
of scattered photons on the gate performance. A pair 
of narrow-waist (< 2 /im) adjacent laser beams in the 
standing- wave configuration produces the necessary force 
to be applied to each ion. Beside this, because of the 
tight focusing it reduces spontaneous emission and laser 
power. Following the treatment of Ref. [Hj], an estimate 
of the needed laser power to realize the logical gate is 
given by 



P = A 



uj xy AhcK 2 w 2 sin 2 (7/2) 



3r|r? 



(34) 



Here A = lul — is the detuning, that is, the dif- 
ference between the laser and the relevant atomic tran- 
sition frequencies, T is the linewidth of the transition, 
k = 27t/Al is the wave number with the laser wavelength 
Xl = 2ttc/ujl, c is the speed of light, w is the size of the 
beam waist, and 7 is the angle between the k vectors of 
the two laser beams (see also Fig. [T]). Additionally, we 
can estimate the influence of photon scattering on the 
gate fidelity as: F scat = e _7Vphot , where the number of 
scattered photons in the standing-wave configuration is 
given by 



phot 



V2TT 3 e cm 2 w 2 Lui y \r? - ff| 3 
3e 2 A L P 



sin 
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FIG. 3. (Color online). Infidelity vs 
0.7, N = 30, and Pg = 4 x 10 3 " 2 
lo c /(2tt) = 76.08 kHz, v xy = u). 
uj z /(2%) = 53.26 kHz, and v T = w r /(27T 
black (solid) line corresponds to r g /r, 
the red (dashed) line to r g /r r = 10" 
the blue (dashdot) line to r g /r r = 10" 



temperature for q z 



£ s mui c . Parameters: v c 
xv /(2tv) = 5.38 kHz, v z 



32.75 kHz. The 
= 10- 1 (r g = 3 /is), 

2 (r g = 0.3 fxs), and 

3 (r g = 0.03 (is), with 



(35) 



r r = 27r/tj r . The inset (on the left corner - top) provides the 
infidelity for r g /r r = 10" 1 (r g = 0.03 /is) with v c = 7.61 MHz 
as in Ref. dgj", v = 300 MHz, v xy = 537.97 kHz, v z = 5.33 
MHz, and v T = 3.27 MHz. The inset (on the right corner - 
bottom) illustrates, for the same trapping parameters as for 
the former inset but with v — 2.4 MHz, the infidelity for 
t s /t t = 10 (r g = 3/is). Such modulation frequency v lies 
within the gap among radial and E x B phonon modes (see 
Fig. El. 



As the last two formulae show, by adjusting 7 we can 
reduce the required laser power, but at the expenses of 
a larger number of scattered photons, and therefore of a 
worsening of the gate performance. 



IV. MODULATED AND STATE-DEPENDENT 
DIPOLE FORCE 



In order to realize our quantum phase gate, |ei, £2) 
e l9eie2 |ei, £2) with ei ; 2 = 0, 1, we have to engineer the 9j~j 
phases in 8 (sec also Sec. [TT|- It is natural to demand 
that the desired value of 8 is obtained with the smallest 
possible value of the applied force (i.e., laser power) or, 
alternatively, in the shortest possible time. This is equiv- 
alent to maximize 8 by maximizing the effect of each 9kj ■ 
This happens when the phases #01 & n d #10 have the op- 
posite sign with respect to the phases 8qq and 6\\. Such 
condition is met when the applied force to the j-th ion 
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FIG. 4. (Color online). Infidelity vs. temperature for q z = 
0.02, N = 30, and P e = 1.3 x 10 5 ^moj c . Parameters: u c = 
7.61 MHz, v xy = 3.80 MHz, v z = 152.16 kHz, and v r = 1.65 
kHz (see text for more details). 

satisfies the relation 

ff ] = -ff ] . (36) 

Additionally, a necessary condition for the implemen- 
tation of a modulated-carrier quantum phase gate is that 
the mean force acting on each ion (respectively each of 
the modes) has to be zero over r g , that is, we have to 
fulfill Eq. (|31l) . Such a requirement can be accomplished 
by making the modulation time symmetric around the 
center of the envelope of the laser pulse. To this aim, we 
impose the further condition on the force: 

£dtf] k) (t)=0 Vfc = 0,l, (37) 

where t = 2-k/v is one period of the modulation. With 
such a condition we obtain a (fast) sinusoidal modulation 
of the force. Experimentally, this can be achieved, for 
example, with an acousto-optical modulator, which can 
vary the frequency of the laser light very quickly. 

A. Energy shifts 

In table U we provide for some ion species the energy 
splitting between the P\/2 and P 3 / 2 levels in the absence 
of an external magnetic field together with the maximal 
value of magnetic field Bz, under which the (normal) 



Atom/Ion 


AE/h (THz) 


Bz (T) 


B PB (T) 


Be II 


1.239 


7.043 


28.170 


Mg II 


17.249 


98.070 


392.286 


Ca II 


41.985 


238.711 


954.845 


Nal 


3.242 


18.410 


73.651 



TABLE I. Energy splitting among the Pi/ 2 and P3/2 levels 
without external the magnetic field (second column from the 
left). In the third and fourth columns (from the left) the max- 
imal and minimal value of the external magnetic field, which 
fix respectively the upper and lower bound for the Zeeman 
and Paschen-Back regimes, are given. 

Zeeman limit can be applied, and the minimal value Bpp 
above which we enter in the Paschen-Back regime. As we 
can gather from the table, the higher the atomic number 
of the ion (or neutral atom) is, the larger AE and the 
limits Bz, Bps- For instance, for the infidelity results we 
showed in the previous section, the corresponding mag- 
netic field at w c = 7.608 MHz arc: B = 4.5 T and B = 12 
T for Beryllium and Magnesium, respectively. These are 
also the values used in current experiments. Thus, for all 
alkaline-earth-metal atoms the Zeeman regime applies, 
and therefore Hp = ^-gjJ z B z well describes the inter- 
action of an ion with the external magnetic field. Here 
gj is the Lande factor [45[ and the nuclear contribution 
has been neglected (gi ~ 10 -3 ). Besides, since the exter- 
nal magnetic field has a strength of few Tcsla, the ionic 
hypcrfinc structure can be also neglected. Additionally, 
we note that in the Paschen-Back regime the transitions 
from the energy (split) ground state (5-level) to the one 
of the excited levels (P) are identical for both ground 
levels when the ion is illuminated with a laser beam of 
a given polarization and frequency. Consequently, the 
dipole force (see Sec. IIVB[) would be the same for both 
states, and therefore it would not be possible in such a 
regime to have state-dependent forces. Instead, this is 
not the case for the broken degeneracy of the S and P 
levels due to the Zeeman effect (see Fig. [6]). 

Finally, we note that given the selection rule s — s' = 
on the quantum number of the spin operator, optical 
transitions between the two levels of S1/2 are not allowed. 
This fact reduces the possibility of undesired flips among 
the qubit pair, and therefore decoherence and dephas- 
ing mechanisms arc strongly suppressed. Henceforth, we 
shall consider the lower energy level of Si/ 2 as the logical 
state |0), whereas the upper one as |1). 

B. Dipole force and dipole matrix elements 

The dipole force is produced by an intensity gradient 
of the laser beam illuminating the atom, which is far 
detuned from the relevant atomic transition, whose levels 
arc ac-Stark shifted. Such an energy shift creates an 
additional potential for the particle. For a two-level atom 
and in the large detuning limit A 3> |0|, the dipole force 
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FIG. 5. (Color online). Modes for the parameters considered in Fig. [3] (top row, oi~ = 0.70, Pe = 4 x 10 3 l^muj c ) and Fig. [4] 
(bottom row, a 2 = 0.02, P e = 1.3 x 10 s f s muj c ). 
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+2/3fi B B 



-2/,fy B B 



P 



3/2 



-2/ibB 



+l/3,i B B 



-1/3MB-S 



Pl/2 



whereas the laser is assumed to be a classical light field. 
Here d a b represents the matrix element of the dipole mo- 
ment operator for the transition |o) = \j = 1/2; rrij) — > 
\b) = \j = 1/2, 3/2; rrij) for a given polarization of the 

electric field £ with strength So, and %(f, t) is the spatial 
and temporal pulse shape. 

The bare (unshifted) dctunings are then defined as: 
= w Ll - w Dl and 5u 2 = wl 2 - w D2 , where w Ll 
and Wi 2 are the laser frequencies. In addition, in order 
to reduce the probability of unwanted photon scattering 
processes, we require that 



-V-bB 



1Mb B\ « \S Ul | « AE |/i B B\ « |i Da | « AE. (40) 

In tabic im we provide the expressions of the state depen- 

+!i b b dent forces for all polarizations of the laser field. Since 

"1/2 j = 1/2 for all relevant transitions, hereafter, for the 

sake of simplicity, we shall denote the reduced matrix el- 

cment by A^Di = M-i/2 1/2 an( l -^d 2 = -^1/2.3/27 where 

m i Mjj> := (j,irij || ef || j^m'j). 



FIG. 6. (Color online). Energy shifts due to the Zeeman 
effect vs. the rrij quantum number. The distances among the 
energy levels are not in scale. 

on the lower energy level reads 

f = -±Vm,f)\ 2 , (38) 

where the Rabi frequency Q(t, f) on the atomic transition 
is given by 

hfl = -d ab ■ £(r, t) = -\d ab \£ Q X (r, t), (39) 



C. Different laser configurations 

As it is evident from the table [Til using only a laser 
pulse either on the Di transition line or on the D2 one it 
is not possible to fulfill the condition (pR)]). To this aim 
we need a second laser pulse (see Fig. [7] on the left) with 
a different detuning from the first pulse. In this section 
we discuss several combinations of the laser polarization 
in order to satisfy both ([55]) and (|37|) . 

1. Pulses with the same polarization 

In this case we first generate two laser pulses with dif- 
ferent frequencies but with the same <j~~ polarization, as 
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Polarization 


I 


J (Mi) 


f\0) /--p, \ 
./ ( U 2j 











Mn 2 £ 2 Vx 2 (r,t) 


Wd/o V X 2 (r,i) 




2h(3d Dl +4ix B B/h) 


4«(5 D2 +^ B B/R) 


4^(35 D2 +5At B B/ri) 


7T 




M Dl £ 2 Vx 2 (r,t) 


A4 D2 £ 2 V X 2 (r,t) 


Mn 2 £ 2 Vx 2 (r,t) 


4h(2^ B B/h-38- Dl ) 


4h(38 Dl +2iJ, B B/Fi) 


2S(^ B -B/fi,-3<5D 2 ) 


2h(35o 2 +n B B/h) 


a + 


2h(4 l i B B/h-3S r>1 ) 





WD 2 £ ! Vx ! Ki) 
4ft(5/n B B/fi-35D 2 ) 


Mn 2 £^X 2 (r,t) 
4h(n B B/h-8 D2 ) 



TABLE II. Dipole forces for all polarizations of the laser fields and internal (logical) states. 






FIG. 7. (Color online). Sketch of the pulse sequence in or- 
der to design the necessary dipole forces to implement the 
modulated-carrier phase gate. The upper lines are red de- 
tuned, whereas the lower lines are blue detuned. In the figures 
the time is in arbitrary units. 



FIG. 8. (Color online). Modulated-carrier gate with the same 
polarization of the laser fields: <r~-polarized (left), and de- 
polarized (right). The distances among the energy levels are 
not in scale. 



in Fig. which corresponds to the first sequence of 
pulses (on the left) in Fig. [Jj Such a configuration of 
lasers yields the following state-dependent forces: 



/ID 

;io> 



(41) 



4h(5n 2 + MB^/ft) ' 

To simplify the notation we make the following replace- 
ments: = -^D i ^oi^ 7 X 2 (^ t) an d B = fisB/h, where 
£oi refers to the electric field strength of either the Di 
(i = 1) or D 2 (i = 2) line. Thus, in order to fulfill p6]> . 
we have to solve the equation 



Xr 



2X r 



(5 D2 + B 3(5 D2 + 5B 3<5 Dl + AB 



(42) 



This can be resolved for both the intensities or the de- 
tunings, so one of them can be considered as a given pa- 
rameter. For the cr + -polarization we obtain an analogue 
equation: 



A'r 



Xr 



2Xr 



5 D2 - B 3fo 2 - 5B 3<5 Dl - 46 



= 0. 



(43) 



As an example, we solve equation (|43p. for instance, for 
the intensities, and be obtain 



(44) 



* Dl = (4g-3,5 Dl )(2(5 D2 - 3g) 
Xr, 2 (Sn 2 -B)(3S D3 -5B) ' 

Such a solution, however, fulfills only the condition (|36j) 
but not the one given by Eq. (f37|) . To this aim we need 
an additional two-pulse sequence, as showed in Fig. [7] on 
the right. Such two pulses can have different strengths of 
intensities and detunings, but they must have the same 
spatial and temporal profile x{r,t) of the first sequence. 
Again, we get, if we solve with respect to the intensities, 
a solution like the one given in Eq. (|44[) . which in gen- 
eral will be different from the solution (j44|) for the first 
sequence of pulses. With such solution we can then easily 
satisfy also the mean zero force condition (|3"T)) by adjust- 
ing the ratios of either the intensities or the detunings. 

In Fig. [9] we display a simple example that shows how 
to achieve the necessary laser pulse sequence. We mod- 
ulate the intensities of the blue (&) and red (r) detuned 
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FIG. 9. (Color online). Designed laser modulation forces in 
order to fulfill both (|37p and (|36[) . The blue detuned laser in- 
tensity with a superimposed square wave signal (top-left) and 
similarly for the red detuned one (top-right). The resulting 
signals are depicted on the bottom-left, whereas switching of 
polarization is given on the bottom-right. 




|0) ■ 



FIG. 10. (Color online). Modulated-carrier gate with differ- 
ent polarization of the laser fields. The distances among the 
energy levels are not in scale. 



4h(36 D2 + 5(j, B B/h) ' 



/ 



|o) 



2h(4fi B B/h - 3(5 Dl ) Ah(S D . 2 + hbB/H) 



, (45) 



whereas for the inverted polarization sequence (Fig. [10] 
on the right) we have 



laser signals like I h ^ r \t) = Jg sin 2 (z^i). The sequence 
starts (t = 0) with both lasers with (repolarization and 
an intensity ratio R + = Iq/Iq given by Eq. (pHf. Then, 
at time t = n/u, the polarization of the two laser fields 
is changed to o~ with another intensity ratio given by 
i?_ = lH II. The ratio i?_ will differ fr om R+, since 
in general the dipole moments are different for the two 
polarizations. Thus, by changing the polarization at 
each minimum of the laser intensity and by choosing the 
proper ratio R± we are able to fulfill the condition ([37)) . 
In order to satisfy the condition (|5rj|) we have to design 
furthermore the ratio of the two successive pulses. This 
is done by multiplying the two intensities I b ( r >(t) with 
square wave signals which are displayed in Fig. [S] (top) 
by the black lines. The resulting pulses arc showed in 
Fig. [5] (bottom-left) whose polarization state is depicted 
on the right lower corner. This procedure ensures that 
both the conditions (f3"T|) and are fulfilled. 



2. Pulses with the different polarization 

The situation in which the laser beams have different 
polarization is depicted in Fig. [TO] If we illuminate the 
ion with red detuned and a+ polarized light and with a 
blued detuned and <r_ beam (Fig. [TU]on the left) we have 
the following state-dependent forces 



= Mp 2 £ 2 Vx 2 (f,t) _ A4 Dl g 2 Vx 2 (r,t) 
1 4H((j, B B/h - 5 D2 ) 2h(3S Dl + 4 f i B B/h) ' 
f\o) = M B2 £^x 2 (r,t) 
J Ah(5fi B B/h-3SB 2 )' 

Both the schemes with the same and with different po- 
larization have the drawback that the transition to the 
excited level P3/2 couples both the ground states of S1/2, 
and therefore producing an additional force that has to be 
compensated with another laser beam. Apart from the 
technical difficulty of putting another laser beam, such 
a beam would also enhance the probability of promot- 
ing an ion to an excited level of Such excitation 
would cause an additional error during the course of the 
gate because of spontaneous emission. Indeed, the ion 
could decay either in the other qubit state or even worst, 
such as for the D-levels in calcium, in another metastable 
state, which would be useless for the purposes of QIP. 

Given that, in order to avoid such scenario, we can 
make still use of the scheme illustrated in Fig. [TUl but 
by avoiding the coupling to the P3/2 manifold, as it is 
showed in Fig. 1111 Here, however, we couple the ground 
state 5i/2 to only the manifold P\/2- The pulse sequences 
are then the same as previously described for the other 
scheme. The detuning from the -P1/2 manifold, however, 
has to be carefully chosen, that is, it has to be much 
smaller than the energy difference among the P\n and 
-P3/2 levels and much larger than fj, B B. Hence, such vari- 
ant works well for sufficiently small magnetic fields. 

In principle there are other possible arrangements ei- 
ther by keeping the frequencies of the laser beams con- 
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FIG. 11. (Color online). Variant of the modulated-carrier 
gate with different polarization of the laser fields. The dis- 
tances among the energy levels are not in scale. 



stant or by keeping constant their intensities. Such com- 
binations rely also on the technical feasibility in an exper- 
imental setup. An important requirement for the design 
of such state-dependent forces is, however, the switch of 
the field polarization, which has to be very fast in order 
to fulfill the condition set by the Eq. (|37[) . as it is also 
shown in the example of Fig. |Hl[d). This can be exper- 
imentally accomplished with Pockcls cells, which can be 
used to manipulate the polarization and the phase of the 
laser. 



V. CONCLUSIONS 

In this work we have analyzed in detail the implemen- 
tation of the modulated-carrier gate presented for the 
first time in Ref. (30| . Firstly, we presented the underly- 
ing idea of the modulated-carrier gate and we provided 
details of the calculations that were only briefly men- 
tioned in Ref. |30j. In that analysis the frame of refer- 
ence rotates at the same frequency of the crystal rotation, 
whose frequency was set to uj c = 2cj r . Within this setting 
the minimal coupling term in the many-body Hamilto- 
nian vanishes. Such approach allows a straightforward 
canonical quantization of the many-body Hamiltonian, 
which reduces to a sum of 37V independent harmonic os- 
cillators. Even though this situation greatly simplifies 
the numerical analysis it does not permit to fulfill the 
condition lo t t s / (2tt) <C 1, which would avoid the uti- 
lization of a co-rotating laser and therefore simplifying 
the experimental realization of the proposed quantum 
hardware. We thus have analyzed the situation in which 
lu c =/= 2uj r . Within this scenario it is no longer possible 
to remove the minimal coupling term in the Hamilto- 
nian of the Coulomb crystal. Nevertheless, by utilizing 
the Williamson theorem for positive definite matrices, we 
were able to diagonalize the classical many-body Hamil- 
tonian, whose normal modes are a combination of both 
the position and momentum variables. As a consequence, 
we were able to perform the canonical quantization. The 
resulting (quantized) Hamiltonian is again given by a 



sum of independent harmonic oscillators. In this new 
situation, however, the matter-field interaction, respon- 
sible of the push on the ion, depends on both conjugate 
"position" and "momentum" operators. We proceeded 
further on by analyzing the performance of the quantum 
phase gate and we showed its robustness for a wide range 
of experimentally accessible temperatures. Importantly, 
we were able to demonstrate that such robustness is also 
displayed for a wide range of ratios T g /r r , therefore allow- 
ing to reduce up to three orders of magnitude the gat e 
operation time compared to the previous analysis [301 ] . 
The drawback is that one has to enhance the modulation 
frequency v up to hundreds of MHz in order to speed up 
the gate operation. We found, however, that by reducing 
the ratio w 2 /u; c , at large values of angular momentum it 
is possible to achieve small rotation frequencies such that 
w r T g /(27r) 1 is fulfilled and high fidelity, for a broad 
range of temperatures, can be obtained with few MHz 
of modulation frequency. This result is quite promising 
since it has been attained with a cyclotron frequency that 
is used in current experiments. 

Finally, we have provided a complete description for 
the design of the necessary forces to be applied on the 
ions in order to accomplish the desired quantum compu- 
tation scheme. To this aim, we have analyzed the ex- 
perimentally relevant region of external magnetic field. 
For all earth-alkali-metal ion species normally used in 
currents experiments the normal Zceman effect provides, 
with good approximation, the right description of the en- 
ergy shifts of the S and P levels. In addition, we have 
also analyzed several possible laser configurations and for 
each one we discussed advantages as well as drawbacks 
and, in some cases, we suggested alternative solutions. 

Further investigations of such a quantum computing 
proposal may rely on further optimization of both the 
force modulation together with a reduced gate operation 
time and its robustness against optimal pulse distortions 
[46| . This can be achieved by means of quantum opti- 
mal control techniques. Beside this, a detailed analysis, 
similar to Ref. @ , in order to characterize and quantify 
all types of errors coming from the quantum dynamics, 
especially due to nonlinearities in the ion-pushing force, 
will be pursued in future investigations. 
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